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Abstract Graphene nanoribbons (GNRs) make up an extremely interesting class of 
materials. On the one hand GNRs share many of the superlative properties of 
graphene, while on the other hand they display an exceptional degree of tunability 
of their optoelectronic properties. The presence or absence of correlated 
low-dimensional magnetism, or of a widely tunable band gap, is determined by the 
boundary conditions imposed by the width, crystallographic symmetry and edge 
structure of the nanoribbons. In combination with additional controllable parame- 
ters like the presence of heteroatoms, tailored strain, or the formation of 
heterostructures, the possibilities to shape the electronic properties of GNRs 
according to our needs are fantastic. However, to really benefit from that tunability 
and harness the opportunities offered by GNRs, atomic precision is strictly required 
in their synthesis. This can be achieved through an on-surface synthesis approach, 
in which one lets appropriately designed precursor molecules to react in a selective 
way that ends up forming GNRs. In this chapter we review the structure-property 
relations inherent to GNRs, the synthesis approach and the ways in which the varied 
properties of the resulting ribbons have been probed, finalizing with selected 
examples of demonstrated GNR applications. 
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1 Introduction 


Graphene has attracted enormous interest since its first experimental realization 
through exfoliation of graphite, mainly because of its many superlative properties 
[1, 2]. By way of example, it is the thinnest, lightest and strongest material known. 
It is also the material with highest thermal conductivity and with highest electron 
mobility. The latter causes graphene to be considered as a potentially revolutionary 
material for future electronic applications [2, 3]. In this respect, however, graphene 
also faces a drawback for its implementation in conventional electronics: the lack of 
a band gap to turn the electron conduction on or off [3]. 

Various strategies have been investigated to overcome this limitation, and one of 
them is through electron confinement in nanostructured graphene. In particular, 
graphene nanoribbons, narrow 1D stripes of graphene with widths in the nanometer 
range, display remarkably varied electronic properties depending on their particular 
structure at the very atomic level [4]. 

A visual framework to understand this is analyzing the aromaticity of some 
representative GNR structures [5, 6]. According to Kekulé, each carbon in poly- 
cyclic aromatic hydrocarbons (PAH) like GNRs or graphene itself has its four 
valence electrons arranged in single, double or triple bonds with electrons from 
neighboring atoms. The structure of a PAH is then the superposition of all possible 
Kekulé bond configurations. Within this picture, the delocalization of six 
m-electrons in a carbon hexagon due to the resonance of two Kekulé configurations 
with alternating single and double bonds is called a Clar sextet, pictured as a circle 
within the corresponding hexagon (Fig. la). According to Clar’s theory [7], the 
most representative and stable structure of a PAH is that with the highest number of 
Clar sextets, which is called the “Clar formula’. Here it is important to keep in 
mind that the bonds sticking out of a Clar sextet are formally single bonds and 
thereby impede two neighboring hexagons to be Clar sextets simultaneously. 

Graphene can be represented with three equivalent Clar formulas, in each of 
which one out of every three hexagons is a Clar sextet, arranged in a (V3 x 73)R30° 
superstructure (Fig. 10-60). Considering a combination of the three possible Clar 


Fig. 1 a Graphical explanation of the Clar sextet as being the superposition of two resonant 
Kekulé configurations. b-d Each of the three equivalent Clar formulas of graphene, with its 
¥3 x ¥3 superstructure highlighted in cyan in (b). e Superposition of all three Clar formulas. 
Dotted circles correspond to Clar sextets within graphene, although not all simultaneously 
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formulas, all hexagons in graphene can be Clar sextets and are fully equivalent 
(Fig. le). In addition, bonds within a Clar sextet are equivalent as well. Since bond 
length alternation (BLA) is a measure of the aromaticity of a system and one of the 
main causes for the opening of band gaps in conjugated organic materials [8, 9], 
graphene is an excellent example of a perfectly aromatic system with absent BLA 
and zero band gap. 

Now we place our focus on GNRs. Because GNRs include many different 
structures, a first classification of GNRs is typically performed depending on their 
edge orientation. They are called zigzag nanoribbons (zZGNR) when the edge runs 
parallel to one of graphene’s lattice vectors, or armchair nanoribbons (aGNR) when 
the edges run along the high symmetry direction thirty degrees off the lattice vectors 
(Fig. 2a). GNRs with edges running along any intermediate orientation are called 
chiral nanoribbons. Their chirality is characterized either by the angle of the edge 
orientation with respect to a graphene lattice direction, or by the edge unit cell 
vector in terms of graphene lattice vectors (Fig. 2a). The key in GNRs is that the 
presence of edges imposes boundary conditions which, among other effects, limit 
the number of possible Clar sextets. 

The width of aGNRs is most commonly given as the number of dimer lines N, 
present across the ribbon (Fig. 2b). In the following we assume the most common 
scenario of purely sp” hybridized carbon atoms within the ribbons, consequently 
with singly hydrogenated edges. Under these conditions, we will see how the 
boundary conditions divide the aGNRs into three families, depending on the 
number of Clar formulas they display [5, 6]. 

For aGNR with N, = 3p, p being an integer, there is only one possible Clar 
formula maximizing the number of Clar sextets. In this formula, all n-electrons 
participate in the Clar sextets, consequently displaying no localized double bonds 
(Fig. 2b). While the n-electrons are delocalized within the Clar sextets, the bonds 
sticking out of them remain single bonds, causing a substantial variation of the 
mean bond length within the ribbons that opens a considerable band gap [5, 6]. For 
aGNRs with N, ಎ 3p + 1 there are two possible Clar formulas, both of which 
include localized double bonds (Fig. 2b). A combination of the two Clar formulas 
still leaves hexagons with and without Clar sextets, causing again notable bond 
length variations within the ribbon. This bond length variation and the electron 
localization in the structure’s double bonds cause this family to have an even larger 
band gap than the 3p family [5, 6]. The scenario is very different with aGNRs of the 
N, = 3p + 2 family. In this case, there are many different Clar formulas, each of 
them featuring two localized double bonds (Fig. 2b). However, a linear combina- 
tion of all the Clar formulas ultimately renders a highly aromatic structure with little 
bond length variations along the middle part of the ribbon. This implies a high 
degree of electron delocalization and consequently a low band gap [5, 6]. 

It is very interesting to analyze zGNRs in this same framework, since they can 
host only a limited amount of Clar sextets even for infinitely long ribbons (Fig. 2c). 
Because sextets can be placed anywhere along the longitudinal axis of the ribbon, 
there are infinite Clar formulas for zGNRs. Their combination ends up being 
equivalent to a superposition of two fully quinoidal structures, that is, structures 
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«Fig. 2 a Classification of GNRs by edge orientation into zigzag, chiral or armchair. The respective 
edge periodicity is given by the unit cell vector in terms of graphene’s lattice vectors a, and ap. 
b Clar sexted distribution on aGNR, classified into three families according to its width in terms of 
dimer lines N,. The family with N, = 3p possesses only one formula, exemplified in 6-aGNRs. 
The family with N, = 3p + 1 possesses two compatible Clar formulas. In each of these formulas 
there is a localized double bond at the edge of the ribbon, as exemplified with the formulas of 
7-aGNRs. On the right, a superposition of both Clar formulas is given, showing all hexagons 
hosting a Clar sextet (although not all simultaneously). The family with N, = 3p + 2 ribbons have 
more than two Clar formulas, namely 3p. Some of them are shown for 8-aGNRs. On the right, a 
superposition of all Clar formulas shows the Clar sextets delocalized along all carbon hexagons. 
c Clar sextet distribution on zGNRs. As exemplified with 4-zGNRs, only a limited number of Clar 
sextets exists even on infinitely long ribbons. Infinite Clar formulas exist, due to their arbitrary 
location along the ribbon, and the superposition is equivalent to the two quinoidal structures 
drawn. Below it is shown how by adding unpaired electrons, Clar sextets can be redistributed 
along the whole ribbon. Some of the Clar formulas are shown, as well as their superposition 
evidencing delocalized Clar sextets along the ribbon and unpaired electrons along the edges. 
d Clar sextet distribution on chiral GNRs, evidencing a limited number of Clar sextets per zigzag 
segment and how addition of unpaired electrons allows creation of additional Clar sextets. One of 
the multiple Clar formulas is shown for each case, including on the right a superposition of all the 
possible Clar formulas that evidences the delocalized Clar sextets along the ribbon and unpaired 
electrons along its edges 


with two double bonds per carbon hexagon (Fig. 2c). However, allowing the 
introduction of unpaired electrons (radicals) into the structure, Clar sextets can be 
distributed anew along the whole ribbon length (Fig. 2c). 

This configuration reduces the overall ribbon energy and readily provides an 
explanation of some of most intriguing properties of ZGNRs. A combination of all 
possible Clar formulas evidences all hexagons as being Clar sextets, implying a 
high degree of aromaticity and electron delocalization along zGNRs, as well as an 
associated low band gap. In addition, because the radicals are always located on the 
edges, combination of the various Clar formulas evidences unpaired electrons fully 
delocalized along both ribbon edges. This is the origin of the edge states in zGNR, 
appearing, as expected from unpaired electrons, at the Fermi level. Finally, since 
unpaired electrons translate to an uncompensated spin, Clar’s formalism also 
implies the spin polarized nature of the edge states in zGNRs. 

Last but not least, chiral GNRs display a combination of alternating zigzag and 
armchair segments along their edges. The lower the chiral angle the longer the 
zigzag segments are, and viceversa. As pictured on the (3,1) chiral nanoribbon 
example in Fig. 2d, there is a limited number of Clar sextets within each zigzag 
segment. However, as occurs for pure zGNRs, addition of unpaired electrons on the 
zigzag edges allows formation of additional Clar sextets (Fig. 2d). The combination 
of the different Clar formulas evidences the high degree of aromaticity and electron 
delocalization in chiral GNRs that explains the low band gap, as well as the creation 
of delocalized and spin polarized edge states along the zigzag segments of the 
ribbon (Fig. 2d). At this point it is important to remind the intricate dependence of 
these edge states with the detailed GNR structure. In infinitely long zGNRs, 
addition of unpaired electrons implies changing from a finite number of sextets to 
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an infinite number of sextets, making the edge state creation always favorable. 
Instead, infinite chiral GNRs readily display infinite Clar sextets (with a limited 
number of them on each zigzag segment, as mentioned above). Addition of 
unpaired electrons simply increases the number of Clar sextets per zigzag segment, 
or in other words the density of Clar sextets per GNR length unit. This has the 
following implications. On the one hand, there is a width threshold below which the 
edge states are not present because for narrow ribbons the energy cost of having 
unpaired electrons may not be compensated by the creation of an insufficient 
number of new Clar sextets. On the other hand, the chiral angle also plays a major 
role. The lower the chiral angle (longer zigzag segments), the more new sextets can 
be created by addition of unpaired electrons, and therefore the lower the width 
threshold for edge state creation. Lastly, the lower the chiral angle the larger is also 
the number of unpaired electrons per edge atom and consequently also the net spin 
moment on each nanoribbon edge. 

A similar analysis based on the Clar sextets and the resulting GNR’s aromaticity 
has been also applied to nanoribbons doped with heteroatoms [10]. Depending on 
the particular atomic species and its bonding site, the Clar sextet distribution and the 
associated aromaticity can change substantially, having a direct and notable impact 
on the GNR’s electronic properties [10]. 

All the above shows, from an intuitive and easy to visualize chemical viewpoint, 
one of the main virtues of GNRs: their amazing variety of electronic properties, 
ranging from large band gap semiconductors to gap-less structures with 
spin-polarized edge states. However, it also remarks the stringent need for atomic 
precision in their synthesis, since minute structural changes at the atomic level can 
cause unproportioned changes in their electronic properties. 


2 Synthesis 


Current digital logic devices require excellent switching capabilities and an on-off 
ratio (Ig,/Iog, corresponding to the current flows under on and off operation con- 
ditions) in the order of 10-107, The latter in turn requires band gaps of 0.4 eV or 
more [3]. If graphene nanoribbons are to be integrated into such structures, their 
width needs to be scaled down to few nanometers [11]. This sets challenging 
boundary conditions to be solved in the synthesis of graphene nanoribbons, com- 
plicated even further by the requisite of atomic precision described in the previous 
section. 

Several groups have readily reported the successful synthesis of GNRs by 
top-down methods and their subsequent characterization. Amongst the different 
strategies we find e.g. the unzipping of carbon nanotubes [12], etching of 
2D-graphene [13], chemical vapor deposition [14], or scanning probe lithography 
[15]. However, in spite of the important insight into the properties of GNRs pro- 
vided in these works, such techniques lack the synthetic reproducibility and/or the 
atomic precision, the latter being of particular importance to be able to aim at 
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rationally chosen GNR structures (e.g. atomically precise width control) and to 
avoid the often dominating disorder effects [16, 17]. 

Instead, Cai and co-workers reported in a seminal work a bottom-up strategy to 
do just that: selectively grow atomically precise nanoribbons with widths in the 
nanometer range [18]. That inspiring work has sparked the interest and research 
efforts in this direction, which basically consists in letting appropriately designed 
precursor molecules to react in a selective way that ends up forming GNRs [4, 19]. 
This approach has been followed in solution [19], as well as supported on solid 
surfaces [4, 19]. Henceforth we will focus on the latter, for which no additional 
functionalization with side-chains is needed to provide solubility to the GNR 
products. In a first step, the reactants are deposited on a surface. Thereafter, they are 
externally activated, triggering polymerization reactions and their subsequent 
transformation into GNRs. Common to all substrate-supported GNRs successfully 
synthesized to date, the reactant deposition was by molecular beam deposition, the 
polymerization based on Ullmann coupling, and the reaction thermally activated. 
Nevertheless, we would like to remark that there may be useful alternatives to each 
of these, as are e.g. electrospray deposition [20], other C-C coupling reactions like 
C-H activation [21] or enediyne cyclizations [22], and photoactivation of the 
reactions [23, 24]. Such alternatives may help reducing the current constraints on 
reactants and substrates. 

A prototypical reaction process is for example that shown in Fig. 3a [25], which 
includes polymerization through Ullmann coupling, but also a second cyclo- 
dehydrogenation step. The whole reaction process has been reviewed extensively [4, 
21, 26] and here we only provide a brief description of the main steps. Ullmann 
coupling itself readily involves various steps [26, 27]. First there is a reactant 
dehalogenation, resulting in surface or adatom-stabilized radical intermediates. This 
process is catalyzed by metallic substrates [27], allowing the dehalogenation to occur 
at temperatures below the temperatures of molecular desorption or of other 
unwanted side-reactions that could compromise the synthetic selectivity and the 
overall polymerization process. In a next step, the radical intermediates diffuse along 
the surface to meet each other and bind covalently. In between the latter two, and 
depending on the type of substrate, the formation of a metastable metal-organic 
intermediate is often observed [21, 26]. Finally, cyclo-dehydrogenation leads to the 
planarization of the polymers, ending up in atomically precise graphene nanoribbons 
with their structure univocally defined by the design of the reactant. Annealing at 
higher temperatures can also cause the coupling of neighboring GNRs by 
inter-ribbon cyclo-dehydrogenation [28-30]. This process, however, results in the 
fusion of a random number of GNRs, lacking selectivity for the formation of precise 
widths and generating substantial disorder. 

Many of the processes described above and their associated energy barriers 
depend notably on the substrate material and surface structure, underlining their key 
role in the overall synthesis. While on the one hand it limits the range of potential 
substrates to be used, on the other hand their appropriate choice can be used as an 
additional parameter to control the reaction and thus the resulting products [31, 32], 
their distribution or alignment [33-36]. 
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Fig. 3 a Schematic representation through wireframe diagrams of an archetypical reaction 
process. The precursors first undergo Ullmann coupling, resulting in non-planar polymeric 
structures due to intramolecular steric hindrance. The circles highlight sections that are elevated 
with respect to the surface. In a second step the planar GNR is formed through cyclodehydro- 
genation. Associated STM images revealing the non-planar and planar nature of polymer and GNR 
are added below the wireframe structures. The use of different precursors allows growing different 
GNRs with varying edge orientation (b—d), width (e-g), edge morphology (h-j) or doping (k—m). 
a Reprinted with permission from [25]. b Reprinted by permission from Macmillan Publishers 
Ltd.: Nature Communications [40], copyright (2013), ¢ Reprinted with permission from [41]. 
Copyright (2016) American Chemical Society. d and i Adapted by permission from Macmillan 
Publishers Ltd.: Nature [18], copyright (2016). e Adapted with permission from [43]. f Reprinted 
with permission from [46]. Copyright (2017) American Chemical Society. g Reprinted with 
permission from [51]. Copyright (2013) American Chemical Society. h Reprinted with permission 
from [45]. j Adapted by permission from Macmillan Publishers Ltd.: Nature [18], copyright 
(2010). k Adapted with permission from [48]. 1 Adapted with permission from [36]. m Adapted 
with permission from [49]. Copyright (2016) American Chemical Society 


Obviously, the particular reactant used also has a notable effect on the various 
chemical processes and their barriers. By way of example, one can compare the 
coupling of 10,10'-dibromo-9,9’-bianthracene [18] and the mono-anthracene equiv- 
alent. The latter lies flat on the surface, and the radical coupling is therefore hindered by 
the steric repulsions between the hydrogens surrounding the carbon radical. Instead, 
the intramolecular steric hindrance between its two anthracene units makes the former 
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a non-planar precursor. As a result, the steric hindrance between two radical inter- 
mediates is reduced, allowing them to get close enough to form the covalent C-C bond 
[4]. 

Another beautiful example is the comparison of the same 10,10’-dibromo-9,9'- 
bianthracene designed to render 7-aGNRs, with 2,2'-dibromo-9,9’-bianthracene, 
designed to render chiral (3,1)-GNRs [25]. For the latter, the Br atoms lie closer to the 
surface, lowering the temperature threshold for dehalogenation. However, even more 
notable is the change in the temperature threshold for the cyclo-dehydrogenation. In this 
case it rather relates to the completely different strain between the two resulting poly- 
mers. With the former, the anthracene units are linked covalently along their short axis 
by a bond that allows free rotational movement with respect to their neighbors. This 
freedom results in alternatively tilted anthracene units along the polymer backbone so as 
to minimize the steric hindrance from opposing H atoms. Instead, with the latter the 
anthracene units are linked covalently to their neighbors both along their long and short 
axes (Fig. 3a). Thus, although the anthracene units still display the same alternative tilt 
to reduce the steric hindrance, the covalent bonds along the long anthracene’s axes limit 
the structure’s rotational freedom, resulting in a substantially strained geometry. 
Sterically induced strain is known to weaken the involved C-H bonds and thereby lower 
the cyclo-dehydrogenation barriers [37—39] in turn explaining the substantially lower 
threshold temperature for 2,2’-dibromo-9,9’-bianthracene [25]. 

Effects like these can be of great importance for potential applications. Focusing 
on the last example, the lower cyclo-dehydrogenation temperature may on the one 
hand allow the use of different substrates that could not stand higher temperatures. 
On the other hand, it may have a strong impact on the resulting products. By way of 
example, on Au(111) it brings the threshold temperatures for Ullmann coupling and 
cyclo-dehydrogenation close to each other. Under this scenario, radical quenching 
by liberated H atoms competes with the radical step growth polymerization, greatly 
lowering the average length of the resultant GNRs [25]. 

However, most importantly, the appropriate design of precursors allows growing 
GNRs with different and tailored structures, all of them with atomic precision. This 
is shown in Fig. 3, displaying selected GNRs with different edge orientations 
(Fig. 30-60) [40-42], different width (Fig. 3e—g) [43-46] different edge structure 
(Fig. 3h—}) [18, 42, 47] and different doping (Fig. 3k—m) [36, 48, 49], all of them 
allowing to tune in a controlled way the electronic properties of GNRs. 


3 Characterization Tools and Associated Insight 


3.1 Electronic Properties Determination 


The electronic structure of GNRs and in particular their band gap (ಗ) size and the 
values of bands’ effective masses (myg and mcg), have been addressed locally with 
scanning tunneling spectroscopy (STS) [34, 40, 42, 43, 45, 50, 51], and with 
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averaging experimental techniques as angle-resolved- inverse- and 
two-photon-photoemission [34, 35, 52-54], high-resolution electron energy loss 
[52, 55] and optical spectroscopy [56]. 

The most widespread method used so far is STS, which measures the local 
density of states (LDOS) as a function of tip position and applied sample bias. It 
allows determining the electronic band dispersion of occupied and empty energy 
levels for the vary same GNR. Precise quantification of bands’ energy level 
alignments and effective masses is given by Fourier-transformed scanning tunneling 
spectroscopy (FT-STS) [57]. The presence of edges in GNRs gives rise to standing 
wave patterns that arise due to scattering of electronic Bloch wave functions against 
such termini. Recording standing waves for different energies along one GNRs 
edge, allows determining the energy versus momentum dispersion of GNRs’ states 
via Fourier transform processing of STS data. Precise bandgap values of 
E, = 2.37 + 0.06 eV and E, = 1.38 + 0.03 eV have been extrapolated from 
valence band maximum and conduction band minimum for 7-aGNR (Fig. 4c-f) 
[501 and 9-aGNRs [45] grown on Au(111), respectively. 
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Fig. 4 Electronic structure of 7-aGNRs. a dI/dV spectrum (red) of a 7-aGNR grown on Au(111). 
The onsets of valence and conduction bands are indicated by the shaded areas. Electronic features 
in the gap are due to the Au surface state (visible in the black spectrum). b ARPES intensity plot 
recorded for 1 ML of 7-aGNRs grown on Au(788) recorded along the ribbon axis. The VB onset is 
found at —0.7 eV below the Fermi level. c STM image and plot consisting on equidistant dI/dV 
spectra of occupied states recorded along the edge of a 7-aGNR. d Line-by-line Fourier transform 
of (€) including a parabolic fit of the VB. The same is done for the CB on a longer ribbon in (e, f). 
a Reprinted with permission from [51]. Copyright (2017) American Chemical Society. 
b Reprinted with permission from [34]. Copyright (2012) American Chemical Society. 
c-f Reprinted figure with permission from [52]. Copyright (2015) by the American Physical 
Society 
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Such precise value could not be simply inferred by single STS spectra (Fig. 4a). 
On one hand, this is due to the fact that the onset of the GNR’s electronic states 
could be ambiguous due to fingerprints of the electronic structure of the substrate. 
On the Au(111) surface, for example, the presence of the surface state at —0.45 eV 
makes it difficult to identify the valence band position in many GNRs. On the other 
hand the signal from GNRs’ valence or conduction bands could be too weak to be 
detected by STS due to the symmetry of the GNRs’ orbitals and the LDOS decay 
above the GNRs’ plane [45, 50]. 

The occupied frontier band of GNRs have been measured in several cases by 
angle resolved photoemission (ARPES), which provides access not only to the 
bands’ onset energy, but to the whole energy versus parallel momentum dispersion 
of occupied GNRs states (as displayed in Fig. 4b) [34]. Being an ensemble aver- 
aging technique, domains of equally oriented GNRs are needed to get k-resolved 
ARPES data. This has been managed by growing the GNRs on vicinal surfaces that 
exhibit a periodic array of steps, which act as templates and drive the uniaxially 
oriented growth of GNRs along the terraces direction. Slight deviation between the 
value of the valence band maximum in flat and vicinal surfaces are found due to the 
different substrate’s work functions. 


3.2 Chemical Structure 


X-ray photoemission spectroscopy (XPS) is the technique of choice to study the 
chemical composition of a material. In the case of GNRs, it has been used to shade 
light on the chemical mechanisms underlying the three basic steps in GNRs’ for- 
mation by on-surface synthesis. After room temperature deposition of precursor 
molecules on the substrate, the shifts of the core levels (as Br3d and Cls) are 
measured as a function of increasing sample temperature in different experiments 
(see Fig. 5 for an example). This allows obtaining invaluable information about the 
whole reaction process including the step sequence, the halogen desorption from 
the surface and the associated threshold temperatures [25, 30, 58]. 

The chemical structure of GNRs has been imaged in UHV with STM and with 
atomic resolution by non-contact atomic force microscopy (nc-AFM). In molecular 
imaging, STM is sensitive to the density of states near the Fermi level, typically 
delocalized over the entire molecule. The correlation of bonding structure and STM 
contrast is thus not straightforward. However, simulations of STM images, reliably 
performed with DFT, are of great value for the interpretation of the STM data and 
finally allow identifying the structure and adsorption geometry of GNRs 
[18, 43, 46]. 

The AFM technique allows, under specific conditions, sensing the short range 
forces between the probe and the measured sample, whereby the signal becomes 
sensitive to the total charge distribution on the adsorbates (which is largest on 
atomic sites and along chemical bonds). Although not strictly required [59, 60], 
such imaging is most commonly performed at low temperature (~5 K), in 
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Fig. 5 Photoemission spectra of the C 1s core levels of 2,2’-dibromo-9,9'-bianthracene deposited 
on (a) Au(111) and (d) Ag(111) held at room temperature and their evolution as a function of 
sample annealing temperature. Similar measurements of the Br 3p core levels are shown in panels 
(b) and (e). Panels (c) and (f) depict Br 3p spectra, together with their associated fits (blue and green 
lines correspond to organic and metal-bound Br components, respectively), of two representative 
temperatures marked with the colored lines in (b) and (e), respectively. The spectra are shifted along 
the intensity axis for better comparison. Reprinted with permission from [25] 


frequency modulation and non-contact modes, by keeping an inert probe (typically 
functionalized with a single CO molecule) close enough to the sample [61]. Doing 
so and tracking the changes in the frequency shift of the oscillating cantilever force 
sensor in a plane above the GNR, it has been possible to precisely resolve the 
bonding structure of pristine and doped GNRs (Fig. 3b—d, k) [40-42, 48], as well as 
defects [40] and atomically sharp junctions [29]. 
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Besides its imaging capabilities, nc-AFM has been used also to proof funda- 
mental properties of GNRs as their structural superlubricity on gold surfaces [62]. 
GNRs with lengths between 5 and 55 nm could be moved laterally by the force 
sensor tip with static friction forces between 2 and 200 pN. Such ultralow forces 
arising in the sliding motion of GNRs pave the way for creating frictionless 
coatings. 


3.3 Detection of Vibrational Modes 


The GNRs vibrational structure predicted by theoretical studies, based on first 
principle methods, is characterized by graphene-like and intrinsic modes [63]. The 
graphene E>,-like (or G) mode is induced by the relative motion of neighboring 
atoms and tends to the value of 1580 cm of the E>, in graphene as the GNRs’ 
width increases. A localized mode is found at 3000 cm" which is the typical 
vibration of the C-H bond. The most prominent peak in the low frequency range 
(200-800 cm") is the radial breathing mode (RBLM) due to the outward motion 
of the nanoribbon’s edge atoms by keeping the central atoms at rest (inset in 
Fig. 6b). This is the most representative vibrational mode of both armchair and 
zig-zag GNRs and since its frequency is roughly proportional to the inverse square 
root of the GNR’s width, it can allow identifying the types of ribbons present on a 
surface. 

Experimentally, the frequencies related to such vibrational modes of GNRs have 
been determined by Raman spectroscopy. Raman spectroscopy measures the 
inelastic scattering of photons by phonons, usually by means of a laser as photon 
source tuned in a spectral range between the near infrared to the near ultraviolet. 
Raman spectra measured on a layer of 7- and 9-aGNRs grown on Au(111) iden- 
tified the different modes predicted by theory, as the RBLM mode at 396 cm" for 
7-aGNRs and 312 cm’! for 9-aGNRs, but also the D peak related to disorder in 
graphene (Fig. 6 a, b) [18, 45]. Raman peaks showed a strong dependence on the 
excitation source wavelength so that the RBLM of 9-aGNRs is depleted by incident 
green laser irradiation (2.33 eV) but enhanced by infrared light (1.58 eV), the 
opposite occurs for 7-aGNRs. According to tight binding calculations coherent 
radial-breathing-like phonons in GNRs are excited for photoexcitation near the 
optical adsorption edge [64]. In fact, in the case of 7-aGNRs, aligned uniaxially on 
Au(788), the size of the optical gap is 2.1 eV as determined by reflectance differ- 
ence spectroscopy measurements (RDS) [56]. Such value is then close to the green 
light excitation while infrared lies within the gap. 

Raman is a powerful technique for the characterization of GRNs under several 
aspects: (1) it does not need ultra-high-vacuum conditions to be operated; (11) it can 
be used to detect GNRs vibrational fingerprints also when they are deposited or 
grown on insulating surfaces; (iii) it allows identifying the orientation of the GNRs 
edges and the GNRs width [65-68]; (iv) it helps to evaluate the quality of a GNRs 
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Fig. 6 a, b Raman characterization of 9-aGNRs grown on Au(111). a Experimental Raman 
spectra taken with two different excitation sources, 785 nm (solid red) and 532 nm (solid green) 
laser lines. Corresponding DFT based simulations done for photon energies of 1.58 eV (dashed 
red) and 2.33 eV (dashed green). b Theoretical dependence of the wavenumber of the RBML with 
the aGNRs width and experimental data for 7-aGNRs and 9-aGNRs. c, d Vibrational HREEL 
spectra corresponding to (c) linear polymer obtained after dehalogenation of DBBA at 470 K on 
Au(111) and (d) 7-aGNRs obtained after cyclodehydrogenation at 670 K. a, b Reprinted with 
permission from [45]. Copyright (2017) American Chemical Society. c, d Reprinted figure with 
permission from [52]. Copyright (2012) by the American Physical Society 


sample by quantifying the defect density via D band [65]. Ideally it will also allow 
to analyze GNRs length distribution. 

The “electron analogue” of Raman spectroscopy, naively, is high-resolution 
electron energy loss spectroscopy (HREELS). This technique uses the inelastic 
scattering of low energy electrons to measure the discrete vibrational energies of 
vibrational modes of adsorbates on a surface. Angle resolved HREELS has been 
used to follow and characterize the thermally activated steps of GNRs formation 
[52, 55]. For GNRs well defined changes in the vibrational spectrum from the 
polymeric phase to the aromatic structure have been measured. The quenching of 
vibrational modes with dipole moments laying within the phenyl rings plane 
demonstrated the planarization of all the phenyl rings after cyclodehydrogenation 
(Fig. 6c, d). 
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4 Tuning the Electronic Properties 


Among the most attractive virtues of graphene nanoribbons is the remarkable 
dependence of their electronic properties on the detailed GNR structure. As a 
consequence, even minimal structural changes brought about in a controlled way 
through the synthesis process, allow tuning the nanoribbon’s electronic properties 
over a wide range. In the following we review some of the different routes that have 
been proposed and explored to control the atomic structure of graphene nanorib- 
bons, as well as their effect on critical electronic properties as for example the band 
gap, band dispersion or energy level alignment. 


4.1 Tuning Through Edge Orientation 


As readily discussed in the introduction in the light of Clar’s theory, the edge 
orientation has a critical effect on the nanoribbon’s electronic properties. On a more 
quantitative basis, we will now discuss some of the details in reference to the 
electronic properties of 2D graphene. Tight-binding calculations have been suc- 
cessfully used to describe with reasonable accuracy the electronic properties of 2D 
and nanostructured graphene. Electrons involved in gaphene’s o-bonds are more 
strongly bound and thus the main interest for optoelectronic processes arises from 
the delocalized m-electrons. Consequently, each atom can be simply characterized 
by a single 2p, orbital. Applying the tight-binding Hamiltonian at the nearest 
neighbor level to the two-dimensional honeycomb structure of graphene with its 
two atoms per unit cell (each forming an equivalent hexagonal sublattice, Fig. 7), 
the following dispersion relation is obtained [69]: 


Ky ky ka k 
E(k) = ay] 3-+2e0s( 5 + = +2cos (“Ss — = +2cos(kya) (1) 


where s = + 1 (—1 for the valence and +1 for the conduction band), t is the hopping 
integral and a is graphene’s lattice constant. Near the 1' point, both valence and 
conduction bands depend quadratically on k, and k,. At the M points there is a 
saddle point in the energy dispersion. The most interesting point is around the K 
points, where the bands show a linear dispersion in so-called Dirac cones that 
display electron-hole symmetry (Fig. 7c) [69, 70]. That is, this first order 
tight-binding calculation already predicts the semimetallic, gapless nature of gra- 
phene, as well as the linear dispersion around the K points in a massless electron 
behavior. 

Applying the same theoretical analysis to GNRs of arbitrary edge orientation, the 
following boundary condition is set: because H atoms along the edges do not 
contribute any m-electron, the GNR’s n-electron wave functions on them vanish. In 
aGNRs, this boundary condition is fulfilled by eigenstates of graphene with 
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Fig. 7 a Graphene lattice and its primitive vectors a; and a. b Reciprocal space of Graphene and 
its reciprocal primitive vectors bj and 0). The Dirac points are located at K+. Dashed lines indicate 
the growth directions of armchair and zig-zag GNRs in the reciprocal space. € Band structure of 
graphene. The close up shows the linear behavior of the bands close to the Dirac points. The slope 
of these cones correspond to the Fermi velocity vr. € Reprinted figure with permission from [70]. 
Copyright (2009) by the American Physical Society 


perpendicular momenta k, that allow standing waves perpendicular to the aGNR 
axis. A detailed derivation is provided in reference [69], resulting in the discrete 
values 

r 20 


= Z r=1,2,3,...N 2 
೫.12 2೩೯ ₹6 (2) 


all 


The associated aGNR’s band structure is thus simply a combination of slices 
across graphene’s band structure at N equidistant kų values. This is shown by way 
of example in Fig. 8 for 9-aGNRs. The frontier bands stem from the cut closest to 


E/ Htl 


Fig. 8 a First Brillouin zone of graphene (black hexagon) and its m-band structure (black/gray 
surfaces) with inset displaying bands near the Dirac cone. Also shown are the first Brillouin zone 
(red rectangle) and cutting lines (red dashed lines) for infinite length 9-AGNRs. b Energy bands of 
empty and occupied states corresponding to cutting lines r = 5, 6,..., 9 in (a). Reprinted with 
permission from [45]. Copyright (2017) American Chemical Society 
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graphene’s Dirac point at K. Depending on the nanoribbon’s width, the properties 
can thus vary from a metallic to a semiconducting behavior. The former occurs with 
aGNRs of the N, = 3p + 2 family, which always present one cut right through the 
K point, while the former applies to the other two subfamilies. At this point it is 
important to note, however, that this tight-binding derivation is based on equivalent 
carbon atoms throughout the nanoribbons. Taking into account that carbon atoms at 
the edges are saturated with hydrogen, such equivalence is lifted and the resulting 
edge distortion actually results in the opening of a band gap also for the 3p + 2 
family [11, 71], as obtained also from more refined ab initio calculations [11, 71] 
and as readily proved experimentally with 5-aGNRs [43]. Further experimental 
examples reporting the electronic properties of different atomically precise aGNRs 
and confirming the great variety of band gaps as a function of their width will be 
discussed in the next section. 

The case of zGNRs is more complicated. As derived in reference [69], the wave 
functions satisfying the boundary conditions have transverse momenta displaying a 
dependence on the parallel momentum such that 


sin(k,N) + gx sink (N+1)]=0, gx 72006 (2). (3) 


That is, most of the bands can still be obtained from slices of graphene’s band 
structure, but those slices are now curved, as displayed by way of example with the 
black solid lines in Fig. 9a for a zGNR with N, = 4 (4-zGNR). It is immediately 
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Fig. 9 a Reciprocal lattice of graphene showing the discrete values of k, for N = 4. The cuts are 
not linear, in contrast with aGNR. Also note the change in orientation of ಓ and k | with respect to 
Fig. 8. The red dotted lines mark the regions where the additional edge states appear that cannot be 
traced back to the band structure of graphene. The reciprocal unit cell of a 4-zGNRs is also shown 
in black dashed lines. b Band structure of a 4-zGNR, revealing the degenerate edge states as k || 
approaches 7/೩. c Band structure of a 50-zGNR revealing the flat bands of the edge states even 
more clearly. d Graphene’s band structure projected into TK direction, showing the absence of the 
flat bands of zGNRs. 7/೩ indicates the boundary of the first 1D Brillouin zone of the GNR (black 
dotted lines in (a)), where the zone-folding technique should be applied. b Reprinted figure with 
permission from [73]. Copyright (1996) by the American Physical Society. € and d Reprinted 
figure with permission from [72]. Copyright (1999) by the American Physical Society 
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obvious that under this tight-binding approximation zGNRs are gapless because one 
of the slices always goes through graphene’s K point. This is also evident from 
Fig. 9b, which displays the band structure associated with these ribbons. It is 
interesting to see in Fig. 9b that valence and conduction band join into a degenerate 
flat band above a critical k || value k. It coincides with the reciprocal space region 
where only N — 1 slices are present instead of N (Fig. 9a), because no k; are found 
as solutions to Eq. (3). In that region, the states giving rise to the flat bands can be 
understood as imaginary solutions to Eq. (3) atk, =Oandk, = 7/೩ as displayed 
with dotted red lines in Fig. 9a. The wave functions associated to those solutions 
are localized at the nanoribbon edges and strongly decay toward the interior of the 
ZGNRs [69, 72, 73]. It is best observed that these edge states cannot be derived 
from slices in graphene’s band structure comparing its projection along the TK 
direction (Fig. 9d) with the band structure of wide z-GNRs (Fig. 9c), which has a 
large density of slices and could thus be expected to be comparable to that of 
graphene. Such comparison immediately makes obvious that the flat bands asso- 
ciated to the edge states are not present in 2D-graphene, but that they are a property 
of zGNRs intrinsic to their zigzag edges. 

The strong density of states at the Fermi level caused by the edge states should 
induce either a lattice distortion via electron-phonon interactions or magnetic 
polarization via electron-electron interactions [69]. In the case of GNRs the latter is 
favored and drives the spin polarization of the edge states. Edge sites along a ZGNR 
are made up by atoms of the same sublattice on one side and of the other sublattice 
on the opposite side. Thus, the edge state being nonzero only on one of the two 
sublattices at each edge, the magnetic moment selectively increases on that sub- 
lattice with a ferromagnetic spin configuration and decays toward the ribbon inte- 
rior. The opposite edge sites belong to the other sublattice and display 
magnetization with opposite spin. The edge states are thus not only spin polarized 
ferromagnetically along each edge, but additionally antiferromagnetically coupled 
across the ribbon [69]. Notably, the electron-electron interactions driving the 
magnetization concurrently open a band gap A° between the edge states, deterring 
ZGNRs from truely being gapless structures. This can be observed in the calculated 
band structure of zGNRs with and without electron-electron interactions in 
Fig. 10b. While tight-binding density of states has only one van Hove singularity 
related to the edge state’s flat bands at E = 0, including the electron-electron 
interaction U term in a mean field Hubbard model results in two pairs of density of 
states peaks split by A° and A! (Fig. 10b). In particular, A° relates to the antifer- 
romagnetic correlation between the two edges, while the larger splitting A' relates 
to the ferromagnetic correlation between the spins of the most strongly localized 
states at ಓ|| = 7/೩ in the same edge. It should be mentioned here that in the 
nearest-neighbor approximation the electron-hole symmetry remains unchanged, 
but it is lifted as a next-nearest-neighbor hopping term is included in the Hubbard 
model [74]. The opening of a bandgap has been confirmed on atomically precise 
6-zGNRs synthesized on Au(111) that were subsequently transferred onto NaCl 
bilayer islands, indeed revealing under such circumstances relatively large splitting 
values of A° = 1.5 eV and A’ = 1.9 eV [42]. Although with somewhat disparate 
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Fig. 10 a Band dispersion obtained from tight binding calculations of GNRs of different 
orientation, from zigzag, through chiral, to armchair. The closer to the armchair direction, the 
lower the density of states of the flat bands at E = 0, until a band gap opens for armchair ribbons. 
Similar calculations for zigzag and chiral ribbons are shown in (0-6) upon addition of 
electron-electron interactions, revealing the opening of a bandgap also for these orientations. 
Reprinted figure with permission from [76]. Copyright (2011) by the American Physical Society 


values, similar conclusions have been drawn also from experiments on ZGNRs with 
worse defined and less perfect edges obtained by differen methods like hydrogen 
etching of 2D-graphene [13] or scanning probe lithography [15]. 

From a structural point of view, chiral ribbons combine zigzag and armchair 
segments. From an electronic point of view the zigzag segments are most deter- 
mining, providing chiral ribbons with low band gaps and spin-polarized edge states. 
However, these properties vary with the chiral angle (0, Fig. 2) in a relatively 
smooth way in between the two limits of zigzag (0 = 0°) and armchair graphene 
nanoribbons (8 = 30°). In the infinite-width limit, the band structures of chiral 
GNRs can be obtained from a continuous rotation of the band structure of graphene 
and the average density of edge states per edge length po at the Fermi level is given 
by [75] 


ae = cos(0+ ಸ). (4) 
30 3 

As discussed above, the large po is responsible for the magnetization of GNRs. 

And as evidenced in Eq. (4), Po is highest for lowest chiral angles and vanishes as 

the armchair angle is approached [74, 76]. This has the important implication that 

the magnetic moment per length unit increases as the chiral angle is reduced. 

Tight-binding calculations of chiral GNRs of finite width as a function of the chiral 
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angle confirm this picture [76]. For GNRs with edge orientations close to aGNRs 
the edge states associated to the flat bands at the Fermi level are almost completely 
suppressed (Fig. 10a). Instead, in ribbons with edge orientation close to a ZGNR the 
flat edge state band extends over the whole 1D Brillouin zone and becomes multiple 
degenerate due to band folding (Fig. 10a) [76]. As in zGNRs, electron-electron 
interactions split the degeneracy of the flat bands also in chiral GNRs (Fig. 10c, d). 
As discussed earlier, each pair of peaks arising in the density of states has a different 
nature, whereby A° and A! relate to the magnetic coupling across the nanoribbon 
and along its edge, respectively. As can thus be intuitively expected, for a given 
GNR width A’ hardly changes with the chiral angle. Instead, A’ decreases in a 
similar way to the net magnetic moment on each edge, as the chiral angle increases 
and approaches the armchair orientation [76]. 

Lastly, it should be remarked that the density of edge states as a function of 
chiral angle depends at the same time on the nanoribbon width. Beyond the visu- 
alization through the Clar sextet theory outlined in the introduction, this effect has 
been computed calculating the number of states at E ಇ 0, normalized to the total 
number of states in the GNR [73, 77]. Calculations reveal that there is a critical 
width below which there is no notable edge state density. Importantly, that critical 
width decreases with decreasing chiral angle (as the edge orientation approaches 
that of zGNRs) [73, 77]. In other words, for a given width there is a minimum 
number of zigzag sites (in relation to armchair sites) along the ribbon’s edges to 
host a notable edge state density. The required ratio of zigzag to armchair sites 
becomes higher for narrower GNRs. 

It is important to observe at this point that the predictions described above for the 
electronic structure of chiral GNRs still lack experimental proof on atomically 
precise chiral GNRs. To date, only the precise synthesis of narrow (3,1)-GNRs has 
been published [25, 31, 32] and only few details of its electronic properties are yet 
known [65]. There are few other works reporting the electronic structure of chiral 
graphene nanoribbons with worse defined edges obtained from unzipped carbon 
nanotubes [12] or chemical vapor deposition [14]. In either case the opening of 
band gaps through the magnetism-driven splitting of the edge states has been 
confirmed by scanning tunneling spectroscopy measurements. 


4.2 Tuning Through Width Control 


Besides the edge orientation, the electronic structure of graphene nanoribbons also 
displays a marked dependence on their width. It is most straightforward to 
understand this dependence for armchair graphene nanoribbons, which is further- 
more the family with the largest number (albeit still very limited) of successfully 
synthesized examples on which to experimentally corroborate such behavior. We 
remind that the band structure of N-aGNRs can be pictured as a combination of 
1D-cuts across graphene’s band structure at N equidistant kų values, namely at 
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k =y 2 r being an integer between | and N. Keeping in mind that graphene’s 
vanishing band gap and the massless electron behavior occur right at the K-point, 
the closer those 1D-cuts get to the K-point, the lower the aGNR’s band gap and the 
effective masses of its frontier bands will be. That is, some of the aGNRs’ attributes 
can be readily understood from simple geometric considerations: along the 
1D-projection, K and K’ are located at ky = 2೫ and kj = ಸ್‌ (Fig. 8a). It is the 3 
in the denominator that accounts for the presence of three aGNR families, 
depending on whether their width in terms of dimer lines N is equal to 3p, 3p + 1 or 
3p + 2 (p being an integer). There is always a cut going directly through the 
K-point for the 3p + 2 family. This family thus has the lowest band gap and 
effective mass, these values being higher on either of the other two families for 
which no cuts across the K-point are present. However, the wider the GNR, the 
smaller the spacing between the cuts and the closer they get to the K-point. This 
determines the inverse proportionality of band gap and effective mass with 
nanoribbon width within each of the families, although if the width increases atom 
by atom it would imply changes between families that could result in an increased 
band gap and effective mass as well. 

It is worth remarking here that, under the assumption of equivalent carbon atoms 
throughout the GNR, tight binding calculations predict the 3p and 3p + 1 families 
to follow a very similar trend. However, this assumption is inaccurate the closer the 
atoms are to the edge, most evidently for the hydrogen saturated carbon atoms. 
Taking into consideration bond length distortions and the associated changes in the 
hopping constant, important changes appear in the resulting band gap calculations, 
as are a Clear splitting of the 3p and 3p + 1 bandgap trends, as well as the opening 
of a bandgap for the 3p + 2 family. These findings coincide with more refined 
ab initio calculations revealing the band gaps of the three different families to order 
according to E, (3p + 1) > Eg (3p) > Eg (3p + 2) [11, 71]. 

Besides tight-binding [11, 69, 71], several other methods have been used to 
quantitatively describe the width-dependent band gap in aGNRs, like for example 
extended Hiickel theory (EHT) [78], density functional theory (DFT) [11, 79], or 
the GW approximation [79]. At this point, it is important to remind that, to date, all 
the experimentally reported band gaps of atomically precise GNRs are from 
surface-supported ribbons. While less important on weakly interacting substrates, if 
more reactive surfaces are used or reactive functional groups are added to the GNR 
structure, strong GNR/substrate hybridizations may occur. Under such circum- 
stances the band gap will be substantially affected and the GNR’s electronic 
properties may bear little resemblance with those of free-standing GNRs. In addi- 
tion, the band gaps of GNRs have been most commonly characterized by ionizing 
techniques such as STS [34, 43, 45, 46] or a photoemission/inverse photoemission 
combination [35]. This implies on the one hand that it is not ground states but 
quasiparticle energies which are probed, and on the other hand that the surrounding 
dielectric medium (both the GNR itself as well as the substrate) can have a sub- 
stantial effect on those state’s energies through polarization-induced screening. 
Among the various theoretical methods mentioned above, the intrinsically present 
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Fig. 11 Calculated bandgap energies for isolated and Au-supported aGNRs of different width for 
each of the three sub-families (3p, 3p + 1, 3p + 2). Calculations are done with the GW method, 
further adding substrate screening for the supported ribbons. Superimposed we find the currently 
available experimental values of Au-supported aGNRs obtained by scanning tunneling 
spectroscopy, providing good agreement. Adapted with permission from [80]. Copyright (2016) 
American Chemical Society 


electron-electron interactions are best taken into account in the GW approximation. 
It is, however, computationally very intensive and has thus only been applied for 
free-standing GNRs. As a result, calculated GW band gaps still overestimate the 
experimental ones [34]. Best results have been obtained when combining the GW 
approximation with a semiclassical image-charge model to account for substrate 
screening [80]. The screening-derived band gap renormalization varies with the 
substrate and the GNR, being in the 1 eV range for aGNRs on Au (Fig. 11) and 
slightly lower for aGNRs on a NaC] bilayer on Au [80]. Its combination with the 
calculated GW band gaps ultimately provides a remarkably good agreement with 
currently available experimental values for 5 [43], 7 [50], 9 [45] and 13-aGNRS 
[46] (Fig. 11). 

The effective mass is another important figure whose relevance stems from its 
inverse relation to the charge carrier mobility [81], a critical parameter for potential 
applications. However, only few studies have characterized effective masses in 
GNRs. From the experimental side, there are measurements for 7-aGNRs and 
9-aGNRs. For the valence band of the former, very disparate values have been 
reported, ranging from ~0.21 me (measured by ARPES [34, 82]) to 0.41 + 0.08 
Me (measured by STS [50]) and even 1.37 m, (measured by two photon photoe- 
mission [53]). For the valence band of the latter, better matching values of 
0.09 + 0.02 m, and 0.12 + 0.03 m, have been reported from ARPES and STS 
measurements, respectively [45]. Interestingly, there have been different works 
proposing analytical relations between the effective mass and other GNR parame- 
ters. For example, Arora et al. proposed the effective mass to be linearly dependent 
on the GNR band gap in a similar way as that claimed for carbon nanotubes: 
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Table 1 Summary of the calculated and experimentally obtained values for the effective mass of 
aGNR’s valence bands given in electron mass (me) units 


Arora et al. [83] Raza et al. Experiment 
[78] 
Width From From ARPES STS 2PPE 
(N) calc. Es exp. E, 
7 0.154 0.21 0.21 0.21 [34] 0.41 + 0.08 1.37 
0.23 [82] [50] [53] 
9 0.074 0.12 0.09 0.09 + 0.02 0.12 + 0.03 — 
[45] [45] 


Experimental E, values for the m* calculation with Arora’s relation have been taken from Refs. 
[45, 50] and for 7-aGNR and 9-aGNR, respectively. 


m* = (E,/11.37) me, with E, given in eV [83]. In turn, Raza et al. related it to the 
aGNR width W through m* ಎ (0.091/W) m, for the 3p family, m* = (0.160/W) me 
for the 3p + 1 family and m* = (0.005/W) m, for the 3p — 1 family (W being the 
distance, in nanometers, between the C atoms on either edge of the ribbon, based on 
a C-C bond length of 1.44 A) [78]. In the following we provide Table 1 com- 
paring the experimentally reported effective masses to the values predicted by either 
of the models mentioned above. It can be seen that there is an excellent agreement 
between predictions and experiments, best of all when taking the experimental 
values obtained from ARPES measurements. 

The width dependence of the electronic properties of GNRs does not only affect 
armchair oriented ribbons, but also ribbons with other orientations. As described in 
the previous section, the edge state magnetization in zigzag and chiral GNRs drives 
the opening of a band gap A° and an additionally split resonance A’ related to the 
inter-edge and intra-edge magnetic coupling, respectively (Fig. 10) [69, 74, 76]. 
Thus, as expected from their respective nature, A! hardly varies with the ribbon’s 
width. Instead, A° decreases with growing width, although in the zGNR case in a 
monotonic way and faster than for aGNRs [74, 75] [76]. This qualitative behavior of 
A° and A! is at least what results from calculations on free standing zGNRs. 
However, while the same trend holds for A° also in calculations of gold-supported 
ribbons, in that case A! shows a surprising increase with ribbon width [80]. 
Experimental confirmation on atomically precise zGNRs of different width is still 
missing, since only 6-zGNRs have been synthesized to date with atomic precision 
[42]. Nevertheless, closely related measurements have been performed on the edge 
states localized at the zigzag ends of 7-aGNRs. For increasing 7-aGNR length, 
which increases the distance between their ends (in analogy to a width increase in 
ZGNRs), the gap between edge states hardly changes. Because the minimum sep- 
aration measured in that study is ~3 nm, this observation is in agreement with a 
rapidly converging band gap at even smaller length scales, further supported with 
calculations [84]. There are also closely related measurement on zGNRs obtained 
from scanning probe lithography and hydrogen etching. The latter shows band gaps 
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Fig. 12 a Flatness index na, defined as the ratio between the number of states with E ~ 0 and the 
total number of states in the ribbon, calculated as a function of zGNR width N,. Density of states 
calculated for zGNR widths of N, = 6 (b), N, = 11 (0) and N, = 51 (d), evidencing the decreasing 
significance of the edge states in the overall density of states for increasingly wide ribbons. 
Reprinted figure with permission from [73]. Copyright (1996) by the American Physical Society 


monotonously decreasing with increasing zGNR width until saturating at widths 
around 4 nm, in reasonable agreement with the predictions described above [13]. 
The former shows a smoothly decreasing band gap with increasing width until a 
threshold value (~7 nm) above which the gap suddenly drops to zero, proposed to 
be related to a sharp antiferromagnetic semiconductor to ferromagnetic metal tran- 
sition [15]. Nevertheless, both studies are based on zGNRs with poorly controlled 
and defined structures as compared to the atomically precise nanoribbons that are the 
focus of this chapter, and thus possibly affected by the additional disorder. 

Beyond the gap, also the edge state magnetization, which is one of the most 
coveted electronic properties of ZGNRs, varies with GNR width. In particular, it 
increases as the width grows, rapidly reaching its saturation value at widths around 
N, = 10 [69, 74]. However, it is important to consider also the relative weight of 
those magnetic edge states on the total density of states of the nanoribbon. Initially 
they increase with increasing N,, but after a peak around N, =7 it ends up 
decreasing at a rate approximately proportional to 1/N, (Fig. 12a) [73]. This effect 
shows how the significance of the special edge state disappears as the ribbon grows 
infinitely wide and becomes 2D graphene (Fig. 12b-—d). 

Chiral ribbons also display a dependence of their electronic properties with their 
width. By way of example, and as readily discussed briefly in the previous section, 
the edge state density in chiral GNRs is absent for widths under a critical value 
[73, 77]. That critical width depends on the chiral angle and decreases as the edge 
orientation gets closer to the zigzag direction. Regarding the band gap, it does not 
decrease monotonically with increasing width as in zGNRs, but displays a periodic 
modulation on top of its overall reduction, reminiscent of aGNRs [77, 85]. 
However, it is worth mentioning that this overall band gap reduction is faster in 
chiral ribbons than in aGNRs. Describing the band gap dependence as E, oc 1/N B, B 
increases from 1 for aGNRs, to larger values as the chiral angle is reduced (e.g. 
B = 1.431 for 0 = 23.41°; B = 2.908 for 0 = 16.1°; or B = 3.758 for 0 = 12.73°) 
[77]. In contrast, the overimposed band gap modulation becomes less pronounced 
for lower chiral angles [77, 85]. 
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Lastly we would like to remind that the discussion above has been made con- 
sidering quasi-infinite ribbons on which only the width is changed. However, 
atomically precise GNRs obtained by on-surface synthesis are finite, typically with 
lengths in the range from few nanometers to some tens of nanometers [25, 35]. As 
for other conjugated molecules or polymers, their band gap scales inversely pro- 
portional to their length, saturating towards the “infinite” limit [9, 84]. Thus, 
although for most applications the growth of GNRs with maximized length is 
desirable, the growth of monodisperse GNRs of limited length could be seen as an 
alternative way to tune and control their corresponding band gaps. 


4.3 Tuning Through Doping 


To understand which effect might doping bring to semiconducting GNRs one can 
look to the effect of impurities on semiconductors [86, 87]. Three scenarios might 
be naively distinguished in: (i) the formation of impurity levels in the semicon- 
ductor’s bandgap; (ii) the appearance of an impurity band; (iii) a rigid shift of the 
whole band structure or modification of the bandgap. In the first case, at light 
impurities concentrations (~10'° atoms cm 3), n-type (or donor) impurities give 
rise to energy levels below the semiconductor conduction band while p-type 
impurities yield electron acceptor levels above the valence band. Doping influences 
the position of the Fermi level, that for n-type moves linearly towards the minimum 
of the conduction band with exponential increase of the doping concentration. In 
the second case at moderate impurity concentrations (~10'* atoms cm” for Si), 
neighboring impurities are so close that their wave functions overlap sufficiently to 
give rise to an impurity band. In heavily doped semiconductors (ಇ2 x 10° cm 
i.e. 3% for Si) the ionization energy of the impurity atoms falls to zero, i.e. the 
impurity band merges with the conduction (for donors) or valence band (for 
acceptors) and only one band is formed. In other cases of heavy doping, as in 
disordered alloys, the band gap can be changed, for example linearly narrowed by 
adding Al atoms in Ga,_,Al,As. 

For carbon-based materials two doping schemes have been used: (i) electrical 
doping which does not alter the lattice structure or chemical composition of the 
material, as by adsorption of a gas or a metal; (ii) chemical doping which can be 
achieved by substitution of carbon (C) atoms with heteroatoms like nitrogen (N), 
boron (B), silicon or sulfur (S). The first type of doping has been reported for 
Li-doped 7-aGNRs [82]. Heavy doping of 7-aGNRs with a monolayer of Li lowers 
the valence band by 1.33 eV and brings the conduction band onset below the Fermi 
level, thereby causing a semiconductor-to-metal transition of the ribbons. Highest 
doping levels (reached upon saturation at around two monolayers of Li) reveal a 
transfer of up to 0.05 electrons per carbon atom. Importantly, such doping is 
accompanied by a quasiparticle bandgap renormalization from 2.4 to 2.07 eV, as 
well as by a dramatic increase of the conduction band’s effective mass. 
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The second type of doping has been successfully demonstrated on atomically 
precise GNRs by means of chemically substituted molecular precursors. The 
inclusion of N, S or B atoms in the GNR matrix corresponds to a heavy doping 
level of several at.% but the resulting electronic structure of GNRs cannot be 
directly related to dopant concentration [88]. In fact, extra-electrons may remain in 
molecular orbitals localized around the dopant heteroatom instead of extending into 
the GNR 7 orbitals [10]. The appropriate design of reactants has allowed to pre- 
cisely introduce dopant heteroatoms either on GNRs edges [36, 55, 89, 90] or 
within their backbone [48, 91]. Ultimately it is the precise combination of het- 
eroatomic species and their localization which defines their effect on electronic 
properties like band gap [10] or energy level alignment. 

In the first case two different situations have been realized. The substitution of 
C-H groups with N atoms in the edges of chevron type GNRs leads only to a rigid 
downshift of the conduction and valence band energies by 0.1-0.13 eV per N atom 
per precursor molecule unit, thus making it a stronger electron acceptor (n-doping 
of the ribbon) [55, 89, 90]. The decrease of the position of the valence band by 
increasing the number of N dopants in the precursor molecule has been also 
characterized with ARPES, STS and DFT for 3-aGNRs grown on stepped Au(788) 
as shown in Fig. 13a—c [36]. The electron lone pair of those nitrogen atoms being 
not in conjugation with the GNR m-system does not further affect the density of 
states of the frontier bands, leaving critical parameters like the band gap or the 
band’s effective mass unchanged. 

On the contrary the substitution of (C-H) groups by S atoms in 13-aGNRs’ 
edges allows hybridization of one of the heteroatom’s lone pairs with the delo- 
calized m-bands. Thus, as predicted by ab initio simulation for free standing rib- 
bons, the energy gap decreases 140 meV with respect to undoped ribbons [49]. 
However, dramatic changes on the energy level alignment are not induced by S 
atoms because their electronegativity difference with respect to C is smaller than in 
the case of N atoms. 

Finally, B impurities have been successfully introduced in the 7-aGNRs’ 
backbone, allowing a full conjugation of empty B p, orbitals with the ribbon’s 
m-system [48, 91]. According to calculations, B-doping induces the formation of a 
new acceptor band lying only 0.8 eV above the valence band. As a result, the band 
gap is strongly decreased with respect to that of pristine 7-aGNR (displayed in 
Fig. 13d, e). Experimentally, the acceptor band and its distribution along the 
B-doped 7-aGNR backbone has been measured by STS, showing good agreement 
with calculations (Fig. 13f, g). However, the position of the valence band (and 
consequently the experimentally determined band gap) has not been reported yet. 


4.4 Tuning Through Strain 


The dependence of the electronic structure of GNRs on the atomic structure of the 
ribbon orientation, the ribbon width and doping have been experimentally 
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Fig. 13 a ARPES spectra of pristine and nitrogen doped poly-para-phenylene (PPP) polymers 
taken at kj close to the valence band maximum. Red lines correspond to parabolic fits of the 
GNR’s bands revealing an unchanged effective mass of 0.19 mọ. b Differential conductance 
spectra performed on PPP (black) and N2-PPP (blue), displaying the rigid shift of valence and 
conduction bands with N-doping. c Computed projected densities of states (PDOS) of the different 
polymers, confirming again the rigid shift of valence and conduction bands upon N-doping and 
underlining the low contribution of the N atoms to the p, states forming the frontier bands. 
d Calculated total density of states (DOS) for pristine 7-aGNRs (dotted line), for B-doped 
7-aGNRs (gray) and contribution from B-atoms to the DOS (red) using the GW approximation and 
including Au(111) substrate screening. e Calculated quasiparticle band structure of B-7aGNRs. 
f Calculated local density of states (LDOS) map of states at the conduction band edge 4 A above 
the borylated ribbon. g Differential conductance map of B-7aGNRs taken at 1 V. ೩-೮ Adapted 
with permission from [36]. d-g Adapted with permission from [91]. Copyright (2015) American 
Chemical Society 


demonstrated by on surface synthesized GNRs. Strain engineering in semicon- 
ducting two dimensional materials offers the possibility to control their optical and 
electronic properties [92]. Thus another route to tune the band-gap of GNRs is to 
exploit strain. To understand the effect of strain on GNRs, it is necessary first to 
revise what occurs on strained graphene. Two types of strain are usually considered: 
uniaxial and shear (Fig. 14a, b). When graphene is stretched in one direction, it will 
shrink in the perpendicular one [93]. Tight binding and first principle calculations 
demonstrated that the effect of a uniform strain applied to the atomic lattice of 
graphene is to drive in the reciprocal space the position of the Dirac cone crossing 
(Ep) away from the K (K’) points while maintaining the cone-like energy disper- 
sion. This Ep will follow a path that is different depending on the type of strain 
applied (shear or uniaxial along armachair or zigzag directions). 

For aGNR shear strain modifies only slightly the band structure but uniaxial 
strain has a relevant effect. The one dimensional Brillouin zone of aGNRs is defined 
by electronic states with allowed k values that lie on parallel lines (k = rz/N + 1 
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Fig. 14 a Schematic view of an aGNR under uniaxial and shear strain. b Calculated energy gap of 
aGNRs as a function of an uniaxial strain £4 on the left, for three ribbons belonging to the three 
families N = 3p, 3p + 1 and 3p + 2 and on the right, for different ribbons of the N = 3p type. 
c Schematic illustration of the effects of strain that induce (i) the deformation of the graphene 
Brillouin zone and motion of the K and K’ points, (ii) a change of the distance between k-lines for 
allowed aGNRs electronic states and (iii) a different crossing of k-lines with the K and K’ points. 
d Behavior of the allowed k-lines near the K point for the three aGNRs families. The green arrow 
indicates the shift of the Fermi point under uniaxial strain. Adapted with permission from [93] 


with r = 1, 2, ..., N + 1) (Fig. 140). The three families of aGNRs characterized by 
their different width (N = 3p, 3p + 1, 3p + 2; p being an integer) have different 
conditions for the crossing of their allowed k-lines with the K (K’) points of the 
Brillouin zone (Fig. 14d). When uniaxial strain is applied Ep moves away from K 
(K’) along the strained Brillouin zone in a direction perpendicular to the k-lines 
(Fig. 140). In this way when Ep arrives in the middle of two lines, the energy gap 
will be maximum; when Ep coincides with ೩ k-line, then the gap will close. This 
finding explains why the energy gap of aGNR is modified in a periodic way with a 
zigzag pattern which is different for each family (Fig. 14b). Moreover, within each 


Bottom-Up Fabrication of Atomically Precise Graphene Nanoribbons 141 


family, the wider the ribbon is, the closer the k-lines are. Thus, the maximum 
energy gap value will diminish. 

For the case of zigzag GNRs, strain can slightly modify the band structure only 
if spin-polarization is taken into account. In this case uniaxial strain will increase 
the bandgap but shear strain will not have much influence. 

To conclude, as observed in the case of graphene, GNR growth on substrates 
which are not lattice-matched might induce strain that eventually is released by 
buckling, fringes, ripples or blisters formation. Interface strain induced on GNRs 
synthesized on a surface has not been observed experimentally yet, but it might 
occur by increasing the interaction of the GNR with its support. 


4.5 Tuning Through Heterostructure Formation 


Heterostructure formation has been successfully used with inorganic semiconduc- 
tors to modulate their electronic properties, whereby sophisticated devices have 
been realized [94, 95]. In analogy, the same concept has been proposed to be 
feasible with GNRs [96, 97], for which one may combine segments with different 
properties like those described above (edge orientation, width, doping, strain), or 
also other parameters like edge terminations or edge structures. Experimentally, 
some examples have been readily demonstrated, creating atomically sharp and 
precise junctions by combination of different reactants (Fig. 15) [89, 98]. An 
alternative way to create randomly distributed GNR heterostructures is by the use of 
only one type of reactant that can, however, react in two different ways [99]. 

First to be realized was a type II heterojunction in which pure hydrocarbon 
reactants were combined with doped reactants that included nitrogen heteroatoms 
(Fig. 15b) [89]. As described in Sect. 4.3 and graphically displayed in Fig. 15b, the 
substitution along the GNR edges of C-H by N has little impact on the electronic 
band gap. In turn, nitrogen’s more attractive core potential lowers the onset energies 
of valence and conduction band and turns the N-doped GNRs into n-type semi- 
conductors. On the other hand, due to the particular GNR and substrate combina- 
tion, undoped GNRs on Au are p-doped, displaying their valence band close to the 
Fermi level. Combination of pristine and N-doped segments into the same 
nanoribbon thus creates sharp p-n heterojunctions, with a band offset of around 
0.5 eV for the heterojunction displayed in Fig. 150, Because the band bending 
occurs over a distance in the order of 2 nm, the resulting electric field at the 
interface is extremely high (2 x 10° V/m), making these heterostructures highly 
promising for electronic devices based on p-n junctions [89]. 

Type I heterojunctions have been synthesized combining precursors that lead to 
differently wide GNRs [98]. The resulting structure is thus a width-modulated GNR 
as shown in Fig. 15a, with segments of 7-aGNRs and of 13-aGNRs. The latter has a 
larger band gap than the former, resulting in the straddling gap evolution that 
characterizes type I heterojunctions. Because the 7-aGNR segments serve as energy 
barrier for charge carriers localized along the 13-aGNR segment, the formation of 
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Fig. 15 a Schematic representation of the precursors leading to differently wide GNRs and an 
associated heterostructure, whose staggered bandgap can be classified as a type I heterojunction. 
Below, constant current STM images display the experimental realization of those heterostruc- 
tures. b Schematic representation of the heterostructure resulting from a combination of pure 
hydrocarbon and N-doped precursors. Below the calculated electronic properties are displayed: 
band structure of pure p-GNRs (undoped precursors), of n-GNR (N-doped precursors), and of the 
density of states of the heterojunction displaying the roughly unchanged bandgap across the 
interface and the type II heterojunction energy alignment. a Reprinted with permission from [98]. 
b Reprinted by permission from Macmillan Publishers Ltd.: Nature Nanotechnology [121], 
copyright (2014) 


quantum well states can be seen for short 13-aGNR segments surrounded by the 
narrower 7-aGNR, establishing an ideal framework for potential band gap engi- 
neering [97]. Furthermore, quantum dots could potentially be used as active 
components inside GNRs. 

An alternative way for the creation of quantum dots embedded inside 7-aGNRs 
has been demonstrated by mixing pristine 7-aGNR with a small amount of boron 
doped 7-aGNRs precursors (Fig. 16a) [100]. The pristine regions in these hybrid 
ribbons preserve the electronic structure of 7-aGNR, while the borylated segments 
lack an energy level aligned with the pristine VB. As a result, the VB electrons on 
the pristine segments become confined by the boron atom pairs, since the VB ends 
abruptly over the borylated regions (Fig. 16b). 

The calculated transmission function (Fig. 16c) for these free-standing hybrid 
7-aGNR, that is, the transmission of electrons from one side of the ribbon to the 
other, shows that the boron pairs are very efficient reflectors and that the VB-1 acts 
as a transmission channel even though the VB is confined (Fig. 16d). This actually 
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Fig. 16 a Schematic representation of the hybrid 7-aGNR, consisting in extended pristine regions 
and disperse borylated segments. b dI/dV maps of a pristine region (dashed blue square) enclosed 
between two borylated segments. The maps show an increasing number of modes in the 
conductance with increasing negative bias, fingerprint of the VB confinement. € Band structure of 
a pristine 7-aGNR (left) and transmission function (right) of a 7-aGNR (red) and a hybrid 7-aGNR 
(blue). The sharp peaks of transmission are related to the quantum well levels caused by the VB 
confinement. The step-like increase of transmission above the VB-1 onset reveals that the VB-1 is 
not confined. d Eigenchannel wave functions of the bands at the energy indicated in panel c. 
The VB electrons are strongly scattered, while the VB-1 electrons transmit almost freely. Adapted 
with permission from [100]. Copyright (2017) American Chemical Society 
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implies that boron atoms selectively confine VB electrons while leaving the VB-1 
unaffected. The reason behind this selectivity stems from the symmetry between the 
boron induced states and the 7-aGNR bands [100]. These results highlight that the 
use of substitutional heteroatoms as dopants, and their related heterostructures, goes 
beyond a simple charge doping model and can strongly modify the transport 
properties of GNRs. 


5 Applications 


Graphene nanoribbons have been used in many types of applications as energy 
storage [101, 102], thin films coatings [103-105] or in composite materials [100-- 
108], but thanks to their attractive physical properties it is in electronic and opto- 
electronic devices applications where they arouse most attention. As a result, 
research efforts have pushed towards the integration of GNRs as active components 
in electronic devices like sensors [109, 110], photodetectors [111-113] and field 
effect transistors (FETs) [114-119]. 

A few experiments have already demonstrated the great potential of GNRs for 
technological applications. For example, nanometer-wide GNRs epitaxially grown 
on silicon carbide were found to be single-channel ballistic conductors at 
room-temperature over distances up to 16 micrometers, which is similar in per- 
formance to metallic carbon nanotubes (which are edgeless but for which the 
chirality determines the electronic structure) [120]. For thinner GNRs, the electron 
transport characteristics have been addressed locally in ultra-high-vacuum by STM 
with single ribbon precision [121]. Transport occurs in a tunneling regime and the 
tunneling decay length through a 7-aGNR measured at different bias voltages 
reveals the dependence of conductance with the ribbon electronic states, that is, 
energy level alignment and band gap value (Fig. 17a, b). 

Graphene nanoribbons technology is still in its infancy and a few problems need 
still to be solved to speed up the use of GNRs in current technology [4]. On-surface 
chemistry, as outlined in the previous pages, could overcome the problems of high 
quality fabrication and functionalization of GNRs. In fact, it allows producing 
graphene nanoribbons with atomically precise control of widths and edges, and thus 
with well-defined electronic structure. It allows to dope GNRs selectively and 
reproducibly, and eventually to increase their functionality by adding active com- 
ponents as photoactive, magnetic or switchable elements. For their operation in 
devices, as in gated multiterminal devices, GNRs have to be transferred to insulators. 
The direct growth on insulating substrates remains a challenge. Promising strategies 
may include the synthesis through photoactivation or through new reactions beyond 
the Ullmann coupling—cyclodehydrogenation combination. Nevertheless, different 
routes can be explored to transfer GNRs from a metal to another material as silicon 
dioxide, as it has been successfully done for two-dimensional materials: exfoliation 
with or without a solution [118, 119], intercalation assisted exfoliation [51], or 
chemical vapor deposition [122]. 
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Fig. 17 a Schematic of the STM pulling experiment (arrow indicates tunnelling current). 
A characteristic current signal during the pulling sequence is shown in the right panel. b Current as 
a function of tip height for different experiments at different bias, evidencing the bias dependence 
of the conductance. ¢ Schematic illustration of the device geometry and scanning electron 
micrograph of the device. d Electrical characterization of a typical device at Vsp = 1 V in both air 
and under vacuum at 77 K. ೩-0 Reprinted by permission from Macmillan Publishers Ltd.: Nature 
Nanotechnology [121], copyright (2012). c, d Reprinted from [118], with the permission of AIP 
Publishing 


The last problem of GNR technology is the final fabrication and processing of 
the device. A few prototypes of field effect transistors have been realized exploiting 
on-surface synthesized GNRs (Fig. 17c, d) [118, 119]. The interest on FETs stems 
on the extreme thinness of GNRs which is expected to enable fabricating FETs with 
very short channels. This should significantly increase the device speed, while 
avoiding the unfavorable short-channel effects of standard CMOS technologies 
[123]. In GNR-FETs the driving current is limited by the contact resistance between 
a single GNR and the metallic contacts caused by the Schottky barrier at the 
interface. The height of the barrier depends linearly with the GNR band gap. In fact, 
FETs with 7-aGNR have shown low driving currents of 1 nA per 1 V drain bias 
[118] while FETs with GNRs displaying lower band gap like 9-aGNR and 
13-aGNR exhibited on-currents up to 100 nA per —1 V bias and on/off ratios of 
10°-10° [119]. To further improve GNR-FETs performance wider ribbons with 
smaller band gaps could be used, or materials for metallic contacts could be tailored 
to the electronic structure of the particular GNR used. In addition, longer ribbons 
(lengths >30 nm) would assure a better contact overlap with source and drain. 
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